function I = trapezoid(f,h)
I=(sum(f(-1:h:1))-0.5*(f(-1)+f(1)))*h;